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1 Introduction 

Konno and Kakuhata are interested in motion of the stretched vortex filaments where the 
filaments are described by the localized induction equation (LIE) LIE is given by 

Xt = X,xXss, (1) 

in which X is the position vector {X,Y,Z), s is the arclength and t is the time. This LIE 
was first obtained by Arms and Hama for the thin vortex filament with the localized induction 
approximation 0. Their original equation is derived as 

^X q X X^ aQ , , 

Under the condition \Xs\ = 1, (O is reduced to (^. is equal to the local stretch of the 

filament Ig defined by [HI 

I -^-\X I 
as 

where I is the length along the filament. If = 1, there is no stretch, if > !> stretch and if 
^5 < 1, then shrink. 

If there is no stretch, one can define the tangent vector t = dX /ds with its magnitude 1 and 
use the Frenet-Serret scheme. Then one introduces the curvature and the torsion, and shows 
that equation is equivalent to the nonlinear Schrodinger equation through the Hashimoto 
transformation |3]. For Ig ^ 1 s does not mean the arclength but the parameter along the 
filament. LIE is an integrable equation. Moreover, we have shown that @ is integrable and has 
a hierarchy which is constructed by the recursion operators [H]. 

In this paper, in order to compare @ with we will try to solve an initial value problem 
with the same initial vortex filaments. The numerical scheme is given in [H]. 

In the next section we show the numerical method for (j2I). In Section 3 numerical results are 
given. The last section presents conclusion. 
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Figure 1. Time evolution (left) and local stretch (right) for the stretched vortex soliton jnj of at 
t = 0, 3, 6 for A = 0.6 and A = 1.5 + i. 

2 Numerical method 

In order to solve we use the following implicit finite difference equation [HI, as an example, 
for the X component as 



At 



1 

I! 



As 



As2 



As 



As2 



(3) 



where 



As2 



The subscripts j and i refer to the time t = jAt and s = iAs, respectively. ^ and Y and Z 
components can be written as the following matrix form 



(4) 



where 



H 



Y. 



D 



Z, 



-Cj,i 

hi 



■J, I 



-Uj^i 



and 



At 



At 



To solve (jlj, we introduce a recurrence formula 



(5) 
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Figure 2. Time evolution (left) and the local stretch (right) for the stretched vortex soliton of 
at t = 0, 3, 6 for A = 0.6 and A = 1.5 + i. 



where 

Ej^i = —(1 — 2L'j^j + Dj^iEj^i^i)~^Dj^i, 

= (1 ~ ^L'j^j + Dj^iEj^i-i) ^{Hj^i — Dj^iFj^i-i), 

for i = 2,3, . . . , N and 

E,, 1 = -(1 - 2D,^i)-^Dj^i, F,- 1 = (1 - 2D,,i)-\Hj^i - Dj^iHj^o) 

Here we take the boundary condition 



H 








As 



(6) 



(7) 



(8) 



The values of all coefficients Q are recursively determined by using Q and ((TJ with Q . Taking 
Xj^N+i = 0, Yj^N^i = and Zj^7v+i = Zj^n + As, we can solve the recurrence formula 

We have taken values of As = 1/8 and At = As^/32. Number of steps j is 12/ At. Numerical 
accuracy is checked against the following aspects, i.e. i) the steady state vortex solution of @ 
with A = 1 propagates without showing significant distortion and ii) the conserved quantity [H] 
J {Is — l)ds conserves within the accuracy 5 % for Q and 10 % for ((^J. 



3 Initial value problem 
3.1 Stretched vortex filament 

As an initial profile of the stretched vortex filament we take modified form of the one vortex 
soliton solution of (0) as 

X = A sin(20 + 5)sech(2G + e), Y = -A cos(20 + 5)sech(2e + e), 

Z = s- ^ tanh(2e + e), (9) 
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Figure 3. Time evolution (left) and local stretch (right) for the loop type of vortex soliton H1U|) of 
at i = 0, 3, 6, 9, 12 for A = 0.5. 



where A is the factor to represent the stretch and we take i = 0. If74 = l, ©is just the exact 
solution without the stretch. $7, and constants 6^ e are given by 



2XrXi 1 , Icq 



2 



ri = \RS — LORt, Q = Xis — uJit, tan(5 = — — 2 -j, e = — -log 



where cq is a constant. The wave number A = A/j + iA/ and the frequency uj = ujh + looi are 
connected by the dispersion relation uj = 2A^. The local stretch is given by 



\2 

sech2(2e + e) - ^ sech^(2e + e) 



li A = 1, there is no stretch. If A 7^ 1, Q shows the local stretch for A> 1 and the local shrink 
for A<1. 

Using LIE, we show the profiles and the local stretch of the vortex filament in Fig. ^ for 
^ = 0.6 with A = 1.5 + i and cq = 1. We observe that, since the filament has shrunk initially, 
then the vortex soliton deforms its profile and moves away with the strong radiation, however, 
the shrunk region remains in the same place as freezing the shrink pp. 

For the case of the equation of motion (jSJ, we illustrate the profiles and the local stretch 
of the filament in Fig. I^Jwith the same parameters such as ^ = 0.6, A = 1.5 + i and cq = 1. 
Comparing with LIE, we observe that the vortex soliton takes a small modification and moves 
fast with the weak radiation. The shrunk region remains in the same place, which is similar to 
the case of LIE. 
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Figure 4. Time evolution (left) and local stretch (right) for the loop type of vortex soliton H1U|) of (|2Jl 
at i = 0, 3, 6, 9, 12 for A = 0.5. 



3.2 Stretched loop type vortex filament 

With the exact solution of as 

X = -A sin 4t sech2s, Y = Acos 4t sech2s, Z = s - tanh2s, (10) 

we consider the initial value problem by taking t = 0. Here A is the factor to represent the 
stretch for the loop type of the vortex filament. The local stretch is given by 

^2 = 1 + 4(^2 _ i)sech22stanh22s. 

In Fig. 121 we show the profiles and the local stretch of the loop type of the vortex filament 
for A = 0.5 at the instants t = 0, 3, 6, 9, 12 by using Since the initial profile is the exact 
solution we observe that the vortex filament rotates in the same place with a constant angular 
velocity and that the shrunk region does not move 

Using the evolution equation we simulate the vortex filament with the same parameter 
as A = 0.5 and show the profiles and the local stretch at the instants f = 0, 3, 6, 9, 12 in 
Fig. 01 From the profiles, we can observe that new propagating vortex soliton is produced by 
accompanied with the radiation. However, from the local stretch, we find that shrunk region with 
double dips remains in the nearly same place, which means that the mode with a constant X, 
a constant Y and two kinks for Z — s excited in the region. 
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4 Conclusion 

We have presented the numerical scheme and considered the initial value problem in order to 
stress the difference of two equations and 0. We have found that the profiles of the vortex 
filaments show quite different behaviour, but for the local stretch, both of the two equations 
give the same results in such a way that the initial stretched region stays in the same place to 
freeze the stretch or the shrink. 
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